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The influence of poor solvent quality on fluid demixing of a model mixture of colloids and nonad- 
sorbing polymers is investigated using density functional theory. The colloidal particles are modelled 
as hard spheres and the polymer coils as effective interpenetrating spheres that have hard interac- 
tions with the colloids. The solvent is modelled as a two-component mixture of a primary solvent, 
regarded as a background theta-solvent for the polymer, and a cosolvent of point particles that 
are excluded from both colloids and polymers. Cosolvent exclusion favors overlap of polymers, 
mimicking the effect of a poor solvent by inducing an effective attraction between polymers. For 
this model, a geometry-based density functional theory is derived and applied to bulk fluid phase 
behavior. With increasing cosolvent concentration (worsening solvent quality), the predicted colloid- 
polymer binodal shifts to lower colloid concentrations, promoting demixing. For sufficiently poor 
solvent, a reentrant demixing transition is predicted at low colloid concentrations. 



I. INTRODUCTION 



Solvents play a crucial role in the thermodynamic 
behavior of macromolecular solutions. Over the past 
half-century, effects of solvent quality on the physi- 
cal properties of polymer solutions have been exten- 
sively studied pi [^]. Polymer-solvent and solvent- 
solvent interactions were first incorporated into the 
classic Flory-Huggins mean-field theory of polymer 
solutions ||. Subsequently, excluded- volume inter- 
actions between polymer segments were identified 
as the key determinants of solvent quality. Polymer 
segments sterically repel one another in a good sol- 
vent, attract in a poor solvent, and behave as though 
ideal (noninteracting) in a theta-solvent. Interac- 
tions between polymer segments strongly influence 
chain conformations and, in turn, phase separation 
and other macroscopic phenomena. 

Compared to solvent effects in pure polymer so- 
lutions, much less is known about the role of sol- 
vent quality in colloid-polymer mixtures. The sim- 
plest and most widely-studied theoretical model 
of colloid-polymer mixtures is the Asakura-Oosawa 
(AO) model |Q, |). This treats the colloids as hard 
spheres and the polymers as effective spheres that 
are mutually noninteracting but have hard interac- 
tions with the colloids. The thermodynamic phase 
diagram of the AO model has been mapped out 
by thermodynamic perturbation theory [^j , free vol- 
ume theory jjj, density functional (DF) theory |^], 
and Monte Carlo simulation By assuming ideal 
polymers, however, the AO model is implicitly lim- 
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ited to theta-solvents. Recently, by incorporating 
polymer-polymer repulsion into the AO model, the 
influence of a good solvent on phase behavior has 
been explored via perturbation theory |Tc| ] and DF 
theory All of these studies assume an effective 
penetrable-sphere model for the polymer coils, which 
is supported by explicit Monte Carlo simulations of 
interacting segmented-chain polymers |l2|, [l3[ B . 
An alternative, more microscopic, theoretical ap- 
proach is the PRISM integral-equation theory jlSj |, 
which models polymers on the segment level. 

The purpose of the present paper is to investigate 
the effect of a poor solvent on the bulk phase be- 
havior of colloid-polymer mixtures. To this end, we 
consider a variation of the AO model that explicitly 
includes the solvent as a distinct component. Specif- 
ically, the solvent is treated as a binary mixture of a 
primary solvent, which alone acts as a theta-solvent 
for the polymer, and a cosolvent, which acts as a 
poor solvent for the polymer. The primary solvent 
is regarded as a homogeneous background that freely 
penetrates the polymer, but is excluded by the col- 
loids. The cosolvent is modelled simply as an ideal 
gas of point-like particles that penetrate neither col- 
loids nor polymers. 

In the absence of colloids, the polymer-cosolvent 
subsystem is the Widom-Rowlinson (WR) model of 
a binary mixture JToj , [l7| , in which particles of un- 
like species interact with hard cores and particles 
of like species are noninteracting. The WR model 
can be shown to be equivalent to a one-component 
system of penetrable spheres that interact via a 
many-body interaction potential, proportional to 
the cosolvent pressure and the volume covered by 
the spheres (with overlapping portions counted only 
once). Hence, in the polymer-cosolvent subsystem, 
the volume occupied by the polymer spheres costs 
interaction energy, inducing an effective attraction 



between polymers reminiscent of that caused by a 
poor solvent. By varying cosolvent concentration, 
the solvent quality can be tuned. Here we investi- 
gate whether and how added hard colloidal spheres 
mix with such effectively interacting polymers. 

In Sec. II, we define more explicitly the model 
colloid-polymer-cosolvent mixture . 



[II 



In Sec. 

develop a general geometry-based DF theory, which 
may be applied to both homogeneous and inhomo- 
geneous states of the model system. The general 
theory provides the foundation for an application to 
bulk phase behavior in Sec. |y|. Readers who are 
inte rest ed only in bulk properties may wish to skip 
Sec. Ill and turn directly to Sec. [Fv|. We finish with 
concluding remarks in Sec. 




II. THE MODEL 

We consider a ternary mixture of colloidal hard 
spheres (species C) of radius Rc, globular polymers 
(species P) of radius Rp, and point-like cosolvent 
particles (species S), as illustrated in Fig. [l| The 
respective number densities are pc( r ), Pp( r ), and 
Ps(r), where r is the spatial coordinate. The pri- 
mary solvent is regarded as a homogeneous back- 
ground for the polymer and is not explicitly in- 
cluded. All particles experience only pairwise in- 
teractions, Vij(r), i,j — C, P, S, where r is the sep- 
aration distance between particle centers. Colloids 
behave as hard spheres: Vcc{ r ) = oo, if r < 2Rq, 
and zero otherwise. Colloids and polymers interact 
as hard bodies via Vcp{r) — oo, if r < Rc + Rp, and 
zero otherwise, and both exclude cosolvent particles: 
Vcs(r) = oo, if r < Rc, Vps(r) = oo, if r < Rp, and 
zero otherwise. The polymers and cosolvent parti- 
cles behave as ideal gases: Vpp(r) = 0, Vss(r) = 0, 
for all r. In essence, this is the AO model with ad- 
ditional point particles that cannot penetrate cither 
colloids or polymers. 

We denote the sphere diameters by <tc = 2Rc 
and up — 2Rp, the bulk packing fractions by r\c = 
AttR c pc /3 and rjp = AirRppp /3, and define a di- 
mensionless solvent bulk density p* s — ps<Jp. The 
polymer-colloid size ratio, q = ap/ac, is regarded 
as a control parameter. 



III. DENSITY FUNCTIONAL THEORY 

We develop a geometry-based DF theory for the 
excess Helmholtz free energy of the model system, 
expressed as an integral over an excess free energy 



FIG. 1: Model ternary mixture of colloidal hard spheres 
of diameter ac, polymer effective spheres of diameter 
ap, and point-like solvent particles. 



density, 

F^ c [pc,PP,Ps] = k B T J d 3 x $ ({<}) , (1) 

where fee is Boltzmann's constant, T is absolute 
temperature, and the (local) reduced excess free en- 
ergy density, is a simple function (not a func- 
tional) of weighted densities, n % v . The weighted den- 
sities are smoothed averages of the possibly highly 
inhomogeneous density profiles, Pi(r), expressed as 
convolutions, 

<(r)= P< (r)*<(r)= / dr' Pi(r>*(r - r'), (2) 



with respect to weight functions, w l u (r), where i = 
C,P,S and v =0,l,2,3,vl,v2,m2. The usual weight 
functions |18, p)) are 



wi(v) = 5(Ri-r), wl{r)=e{R % -r), (3) 



"v2 



(r) 



"m2 



(r) = wi(r) 



rr 



(4) 

where r = |r|, S(r) is the Dirac distribution, 9(r) 
is the step function, and 1 is the identity ma- 
trix. Further linearly dependent weight functions 
are w\[v) = u;«(r)/(47ri2), w^r) = wt 2 (r)/(4^R), 
and u?o( r ) = w\(r)/R. The weight functions for 
v = 3,2, 1,0 represent geometrical measures of the 
particles in terms of volume, surface area, integral 
mean curvature, and Euler characteristic, respec- 
tively Pq. Note that the weight functions differ 
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in tensorial rank: Wq, w\, w\, and w\ are scalars, 
and w* 2 are vectors, and w l m2 is a (traceless) 
matrix. 

The excess free energy density can be expressed 
in the general form 



$ = $c + $cp + ®cs + *CPS, 



(5) 



where the four contributions have forms motivated 
by consideration of the appropriate exact zero- 
dimensional limits. The colloid contribution, $c 
is the same as that for the pure hard-sphere (HS) 
system @ ||: 



$ c = -nfln(l-nf) 
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The colloid-polymer interaction contribution, $pp, 
is the same as in the pure AO case f§ , 



CP 



(7) 



while the colloid-solvent interaction contribu- 
tion 121 is 



-nf ln(l 



^3 C )■ 



(8) 



Finally, in order to model the WR-type interac- 
tion between polymers and cosolvent particles in the 
presence of the colloidal spheres, we assume 



&CPS 



p 

n n 3 

1 — n? 



(9) 



which takes into account the volume excluded to the 
polymer and cosolvent by the colloids. 

It is instructive to compare the current theory to 
geometry-based DF theories previously formulated 
for two related ternary model systems. One starting 
point is a ternary AO model that combines a binary 
HS mixture and one polymer species pif . Letting 
the radius of the smaller HS component go to zero, 
one obtains the cosolvent species. The other start- 
ing point is a recently-introduced model |22| for a 
ternary mixture of colloids, polymers and hard van- 
ishingly thin needles of length L, where the needles 
are ideal amongst themselves but cannot penetrate 
the polymers (hard-core interaction). In the limit 
L — > 0, the needles become identical to the cosol- 
vent particles. We have explicitly checked that the 
DF theories for both systems reduce to the theory 
described above, demonstrating the internal consis- 
tency of the geometry-based approach. 



IV. RESULTS AND DISCUSSION 
A. Bulk Limit 

For bulk fluid phases the density profiles are ho- 
mogeneous: Pi{ r ) = const. In this case, the inte- 
grations in Eq. pft are trivial, and simple expres- 
sions for the weighted densities can be obtained. In- 
serting these expressions into the excess free energy 
density [Eqs. @-@] yields the bulk excess free en- 
ergy in analytic form. The HS contribution, which 
is equal to the Percus-Yevick compressibility (and 
scaled-particle) result, is given as 



37 ?c [3r?c(2 - Vc) ~ 2(1 - Vc ) 2 ln(l - Vc )} 



Vc) 



(10) 

The colloid-polymer contribution is equal to that 
predicted by free volume theory 0, and subse- 
quently rederived by DFT fl : 



rtp/(BirR%) 



{Sqrjc [6(1 - r]c) 



(1 - Vc) 3 

+ 3q(2-r)c-Vc) + 2 q 2 (l + Vc + Vc)] 
- 6(l-?7 C ) 3 ln(l-7,c)}. (11) 

This contribution is linear in the polymer density 
and has a form that arises, as in the original free 
volume theory |Q , from treating the polymers as an 
ideal gas occupying the free volume between the col- 
loids. The colloid-cosolvent contribution is given by 



$CS = -Asln(l - tjc). 



(12) 



This contribution can be similarly interpreted as the 
free energy of an ideal gas in the free volume of the 
colloids. In this case, however, the ideal gas consists 
of point-like cosolvent particles, considerably simpli- 
fying the analytical form of the free volume. In fact, 
by letting q — > in Eq. (11), and identifying species 
P and S, $cp reduces to $cs- The remaining con- 
tribution couples the densities of all three species, 
and is given by 



CPS 



PsVp 



1 



Vc 



(13) 



In the absence of colloids (r/c = 0), this is equiva- 
lent to the mean-field free energy of the WR-model. 
Eq. (13) is a non-trivial generalization thereof to the 
case of non- vanishing r\c- F° r completeness, the re- 
duced ideal-gas free energy is 



$id = 



J2 Pi[MMf)-l], 



(14) 



i=C,P,S 
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where the Aj are (irrelevant) thermal wavelengths of 
species i. This puts us in a position to obtain the 
reduced total free energy density, $ to t = $id + 
of any given fluid state characterized by the bulk 
densities of the three components and the size ratio 

q- 



B. Phase Diagrams 

The conditions for phase coexistence are equal- 
ity of the total pressures, p to t, and of the chemical 
potentials, /Ltj, in the coexisting phases. For phase 
equilibrium between phases I and II, p\ ot = p l ^ at and 
n\ = p] 1 , i = C,P, S, yielding four equations for six 
unknowns (two state-points, each characterized by 
three densities). In our case, a set of analytical ex- 
pressions is obtained from 



Ptot 

k B T 



E • 

i=C,P,S 



d Pl 



and 



Hi = k B T 



dpi 



(15) 



(16) 



the numerical solution of which is straightforward. 

In order to graphically represent the ternary phase 
diagrams, we choose the system reduced densities, 
VCiVPj an d P*s as basic variables. For given q, these 
span a three-dimensional (3d) phase space. Each 
point in this space corresponds to a possible bulk 
state. Two-phase coexistence is indicated by a pair 
of points joined by a straight tie-line. We imag- 
ine controlling the system directly with r\c and r/p, 
but indirectly via coupling to a cosolvent reservoir, 
whose chemical potential, ps, tunes the solvent qual- 
ity. Note that, because the cosolvent is treated as 
an ideal gas, the reservoir's density is simply pro- 
portional to its activity. Thus, the reduced density, 
p*g = exp(/is/fceT), may be equivalently taken as 
a control parameter, which is equal in coexisting 
phases. To make contact with Flory-Huggins the- 
ory, we are implicitly considering here the case in 
which the Flory interaction parameter, Xi falls in 
the range 0.5 < x < lj corresponding to a negative 
excluded- volume parameter, v oc (1 — 2%). 

We initially consider colloids and polymers of 
equal size (ac — crp)- For this case, Fig. || shows 
projections of constant-ptf surfaces onto the three 
sides of the coordinate system, namely the t]c — p%, 
Vc — VP? and r/p — p* s planes, as well as a perspective 
3d view. For reference, the phase diagram without 
cosolvent is shown in Fig. ga. This is identical to 
the common free volume demixing curve of the AO 



model @, |). For p* & T = 0, in which case p* s = 0, 
the r)c — p*s and r/p — p* s planes are inaccessible, 
i.e., all accessible states lie completely within the 
Vc ~ Vp plane. Upon increasing the cosolvent reser- 
voir density to p*g — 0.5, and thus worsening the 
solvent quality, the demixed region grows, as seen 
in Fig. ||b. The critical point shifts towards lower 
7]c and higher r/p, the tie lines become steeper, and 
the area beneath the colloid-polymer binodal in the 
Vc ~ Vp plane (a measure of miscibility) decreases. 
As a physical interpretation of the results, one can 
imagine the polymer spheres as tending to merge 
(overlap) to avoid contact with the solvent. The 
resulting polymer "dimers," "trimers," etc., act as 
larger depleting agents, increasing the range of the 
effective depletion potential between colloids. At the 
same time, the lower effective concentration of de- 
pletants reduces the osmotic pressure and thus the 
depth of the potential. Comparing the phase dia- 
grams for different cosolvent reservoir densities, we 
can conclude that the net effect of merging polymers 
is to increase the integrated strength of the depletion 
potential and thus to promote demixing. 

Eventually, at p* s r = 0.64894, the colloid-polymer 
critical point meets the r/p — p* s plane (where rjc = 
0), as seen in Figs. ^ and (on a larger scale) ||d. 
Polymers and cosolvent here begin to demix already 
in the absence of colloids (the critical point of the 
WR model). For still higher cosolvent reservoir 
densities (beyond the WR critical point), the crit- 
ical point vanishes from the phase diagram and a 
polymer-cosolvent miscibility gap opens up at rjc — 
0. It is tempting to interpret this demixing as ag- 
gregation of the polymer spheres, although it must 
be emphasized that the WR model can only crudely 
describe polymer aggregation. 

Another intriguing prediction is the reentrant 
colloid-polymer mixing evident in Fig. |2|i. For suffi- 
ciently low colloid concentrations and high cosolvent 
reservoir densities (poor solvent), colloids and poly- 
mers initially demix with increasing r/p. Upon in- 
creasing r/p further, miscibility returns over a small 
range before demixing again occurs at higher r\p. 
Such a phenomenon could conceivably result from 
the complex interplay between range and depth of 
the depletion potential arising from solvent-induced 
overlap of polymers. 

For smaller polymer-to-colloid size ratios, the 
above scenario persists. Figure ^ shows qualitatively 
similar results for q = 0.5 and cosolvent reservoir 
densities p* s r = (Fig. |a) and 0.5 (Fig. §3). 
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FIG. 2: Demixing phase diagram of the model ternary 
colloid-polymer-solvent mixture for ac ~ ap and p* s r = 
(a), 0.5 (b), and 0.64894 (c). The latter case is shown 
also on a larger scale (d). 
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FIG. 3: Same as Fig. ^, but for ac = lap and p* s r — 
(a) and 0.5 (b). 



V. CONCLUSIONS 

In summary, we have investigated the bulk fluid 
demixing behavior of model mixtures of colloids and 
nonadsorbing polymers in poor solvents. Our model 
combines the Asakura-Oosawa model of hard-sphere 
colloids plus ideal penetrable-sphere polymers with 
a binary solvent model. The solvent comprises a 
primary theta-solvent and a cosolvent of point par- 
ticles that are excluded from both colloids and poly- 
mers. Cosolvent exclusion energetically favors over- 
lapping configurations of polymers. Although some- 
what idealized, the model exhibits the essential fea- 
ture of solvent-induced effective attraction between 
polymers, mimicking the effect of a poor solvent. 

To study the equilibrium phase behavior of this 
model, we have derived a geometry-based density 
functional theory that combines elements of previous 
theories for the AO and Widom-Rowlinson models. 
Applying the theory to bulk fluid phases, we have 
calculated phase diagrams for cosolvent densities 
spanning a range from theta-solvent to poor solvent. 
With increasing cosolvent concentration (worsening 
solvent quality), the predicted colloid-polymer bin- 
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odal shifts to lower colloid concentrations, desta- 
bilizing the mixed phase. Beyond a threshold co- 
solvent concentration, a reentrant colloid-polymer 
demixing transition is predicted at low colloid con- 
centrations. 

Predictions of the theory could be tested by com- 
parison with simulations of the model. Qualitative 
comparison with experiment also may be possible, 
but would require a relation between the cosolvent 
concentration (as a measure of solvent quality) and 
the Flory interaction parameter. In principle, such 
a relation could be established by calculating the ef- 
fective second virial coefficient of the polymer in the 
polymer-cosolvent subsystem. 

Although here we have approximated the poly- 



mers as mutually noninteracting, their effective at- 
tractions being driven only by cosolvent exclusion, 
future work should include non-ideality between 
polymers, arising fundamentally from excluded- 
volume repulsion between polymer segments. For 
this purpose, a reasonable model is an effective- 
sphere description based on a repulsive, penetrable 
pair interaction (finite at the origin), e.g., of step- 
function or Gaussian shape |l2). The competition 
between such intrinsic repulsion and the solvent- 
induced attraction considered in this work is likely 
to produce rich phase behavior. As a further out- 
look, our approach also could be applied to effects 
of solvent quality on polymer brushes adsorbed onto 
surfaces of colloidal particles. 
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